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Abstract 

We consider mean curvature flow of n-dimensional surface clusters. At (n — 1)- 
dimensional triple junctions an angle condition is required which in the symmet- 
ric case reduces to the well-known 120 degree angle condition. Using a novel 
parametrization of evolving surface clusters and a new existence and regularity 
approach for parabolic equations on surface clusters we show local well-posedness 
by a contraction argument in parabolic Holder spaces. 
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1 Introduction 

Motion by mean curvature for evolving hypersurfaces in M""*"^ is given by 

V = H, 

where V is the normal velocity and H is the mean curvature of the evolving surface. Mean 
curvature flow for closed surfaces is the L^-gradient flow of the area functional and many 
results for this flow have been established over the last 30 years, see e.g. Huisken [T9] . 
Gage and Hamilton [H] , Ecker [9] , Giga [T7] , Mantegazza |2l] and the references therein. 

Less is known for mean curvature flow of surfaces with boundaries. In the simplest 
cases one either prescribes flxed Dirichlet boundary data or one requires that surfaces 
meet a given surface with a 90 degree angle. The last situation can be interpreted as the 
L^-gradient flow of area taking the side constraint into account that the boundary of the 
surface has to lie on a given external surface. A setting where the surface is given as a 
graph was studied by Huisken [20] , who could also analyze the long time behaviour in the 
case where the evolving surface was given as the graph over a flxed domain. Local well- 
posedness for general geometries was shown by Stahl [2S] who was also able to formulate 
a continuation criterion. In addition he showed that surfaces converge asymptotically to 
a half sphere before they vanish. 
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Figure 1: A surface cluster consisting of three hypersurfaces with boundary and one triple 
line on the left; and a surface cluster with four hypersurfaces, where the topology of the 
individual surfaces is not the same for all on the right; taken from [3]. 



Much less is known about the gradient flow dynamics for surface clusters. In this 
case hypersurfaces F^, . . . , F-^ in M""*"^ with boundaries dT^, . . . , dT^ meet at {n — 1)- 
dimensional triple junctions, see e.g. Figure [TJ Here, boundary conditions at the triple 
junction which can be derived variationally have to be described. In what follows we 
briefly discuss how to derive these boundary conditions. We define the weighted surface 
free energy 

N 



^(F) := E / 

-1 



for a given surface cluster F = (F^, . . . , F^) (and constant surface energy densities 7^ > 0, 
i = 1, . . . , N) and consider a given smooth vector field 

C : . 

Then we can define a variation T{e) of F in the direction ( via 

r{e) = {x + eC{x) \xer}. 

A transport theorem now gives 
d 



id'w = - / v'Wdw+ / v^dn 



in—1 



de 

r^(e) r*(e) dV\e) 

where is the normal velocity and is the mean curvature of F*. In addition f * is the 
outer conormal velocity of the surface, i.e. we have w * = (C, i/'), where z/* is the outer unit 
conormal of dV^ (for details we refer to Garcke, Wieland and Depner, Garcke [H]). 
The first variation of T is now given by 



d_ 
de 



..,in.i^un-l 

and hence a suitably weighted L^-gradient flow is given by 



V' = on F* and (1.1) 

3 

E^7V* = at triple junctions. (1.2) 

i=l 
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We remark that the last condition reduces to a 120° angle condition in the case that all 
7j's are equal. 

Local well-posedness for curves in the plane has been shown by Bronsard and Rei- 
tich [6j in a C^"*""'^"^^ setting using parabolic regularity theory and a fixed point argument 
(for a typical solution see Figure [2]) . Kinder lehrer and Liu [21] derived global existence 
of a planar network of grain boundaries driven by curvature close to an equilibrium. 
Mantegazza, Novaga and Tortorelli [25j were able to establish continuation criteria and 
Schniirer et al. |2Z1 and Bellettini and Novaga [1] considered the asymptotic behaviour of 
lens-shaped geometries. We remark that all of these results are restricted to the planar 
case. 




equal area non-equal area 



Figure 2: Mean curvature fiow of a double bubble in the plane, see for results in M^. 

The higher dimensional situation is much more involved as the triple junction now is at 
least one-dimensional and a tangential degree of freedom arises at the triple junction. In 
addition, all mathematical descriptions of the problem result in formulations which lead to 
a free boundary problem. Only recently, Freire [i3\ was able to show local well-posedness 
in the case of graphs. Of course most situations cannot be represented as graphs. We use 
a new parametrization of surface clusters introduced in Depner and Garcke [8] to state 



the problem (1.1), (1.2) as a system of non-local, quasilinear parabolic partial differential 
equations of second order. The PDEs are defined on a surface cluster and are non-trivially 
coupled at the junctions. To simplify the presentation, we will now stick to the situation of 
three surfaces meeting at one common triple junction. But we remark that generalizations 
of our approach to more general surface clusters are possible as long as different triple 
junctions do not meet. Of course this can happen for soap bubble clusters, see Taylor [30] 
and Morgan [26]. In addition we want to remark that in the situation on the left in 
Figure [T] it is in principle possible to use one global parametrization for all three evolving 
hypersurfaces. In this case we would get a system of PDEs on one reference configuration. 
Due to the topological restrictions this is not possible any more in the situation on the 
right in Figure [TJ But since we only use local parametrizations, our method works also in 
this case. 

We hence look for families of evolving hypersurfaces r' (t) C ]R"+^ (i = 1, 2, 3) governed 
by the mean curvature fiow, which is weighted by > (i = 1, 2, 3). These hypersurfaces 
meet at their boundaries as follows 

dT\t) = dV^t) = dV^'it) (=: S(t)), 

which is an (n — l)-dimensional manifold. Also, the angles between hypersurfaces are 
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prescribed. More precisely, we consider 



' V^ = l3'H'' on P(t), te[0,T] (z = 1,2,3), 
Z{r{t),P{t)) = e>^ on tG[0,T], 
((z,j,A;) = (1,2,3), (2,3,1), (3,1,2)), 

[ r{t)\t=o = n (^ = 1,2,3), 



1.3) 



where Tq (i = 1,2,3) are given initial hypersurfaces, which meet at their boundary, i.e. 
OTq = SFq = OTq (=: Sq), and fulfill the angle conditions as above. Here, and if* are 
the normal velocity and mean curvature of r*(t), respectively. 

In (|l.3|), 9^, 9"^ and are given contact angles with < 6'* < vr, which fulfill + + 

sin 6'^ sin 6'^ sin 6'^ 



9^ = 2iT and Young's law 



(1.4) 

^-L 

Let z/*(-,t)(i = 1,2,3) be the outer conormals at dr^{t). Then, introducing the angle 
conditions as in ( 1.3[ ), one can show that (1.4) is equivalent to 



7V^(-,t) +7V'(-,t) +7V3(-,t) = on S(t) 



;i.5) 



which is the condition (1.2) stated above. To choose appropriate normals N'^{-,t) of r*(t), 
we observe that due to the appearance of a triple junction E(t) the six vectors N\-,t), 
z/'(-,t), i = 1,2,3 on E(t) all lie in a two-dimensional space, namely the orthogonal 
complement (To-S(t))"'" of the triple junction. In this two-dimensional space we choose an 
oriented basis and a corresponding counterclockwise rotation R around 90 degree. Then 
we set 

N\-,t) := Riy\-,t) on S(t) 



1.6) 



and extend these normals by continuity to all of r*(t). Then we can write instead of (1.5) 

'y^N\-,t) +-f'^N\-,t) +-f^N^{-,t) = on S(t). 
In the following the angle conditions at the triple line are written as 

{N\-,t),N^{-,t)) = COS9'' 



(1.7) 

on E(t) for k) = (1, 2, 3), (2, 3, 1), and (3, 1, 2). Here and hereafter, ( ■ , ■ ) means the 
inner product in M"+^. 

We are able to show the following result (for a precise formulation of the result we 
refer to Section Isl): 



Main result. 

Let (Fq, Fq, Fq) be a C^"*"" surface cluster with a C^"*"" triple junction curve 7. We assume 
the compatibility conditions 

- (FJ,Fq,Fq) fulfill the angle conditions, 

- -i^(5^Hl + -i'^fi'^Hl + -i^(5^Hl = on the triple line dVl = dTl = dTl 
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Then there exists a local solution of 

= W + angle conditions, 

with initial data (rQ,rQ,rQ). 

The idea of the proof is as follows: First we study the linearized problem around 
a reference configuration with energy methods (this is non-trivial as the system is de- 
fined on a surface cluster). Then we show local (7^"'"°''^"'" 2 -regularity of the solutions to 
the linearized problem. In order to apply classical regularity theory close to the triple 
junction, we parametrize the cluster locally over one fixed reference domain and check 
the Lopatinskii- Shapiro condition for the resulting spatially localized system on the flat 
reference domain directly and for convenience with an energy argument. Finally we use a 
fixed point argument in C*^"'""'^"'" 2^ which is non-trivial as the overall system is non-local. 
In this context ideas of Baconneau and Lunardi |2] are useful. 

We remark that we do not need the initial surfaces Fq to be of class C^"^" as in [2] 
since we linearize around smooth enough reference hypersurfaces, which are close enough 
to Fq in the C^"'""-norm. 

We also remark that the overall problem has a structure similar as free boundary 
problems. This is due to the fact that at the triple junction a motion of the surface 
cluster in conormal direction is necessary. When formulating the evolution on a fixed 
reference configuration, we need to take care of the conormal velocity which results in a 
highly nonlinear nonlocal evolution problem similar as in several free boundary problems, 
see e.g. Escher and Simonett pjj or Baconneau and Lunardi [2]. In our context an 
additional difficulty arises due to the fact that three surfaces who all have a conormal 
velocity meet at the triple junction. The connection to free boundary problems is more 
apparent in the graph case which has been considered by Freire |13j . 



2 PDE formulation 

2.1 Parametrization of surface clusters 

Let us describe F*(t) with the help of functions : F* x [0,T] — )■ M as graphs over some 
fixed compact reference hypersurfaces F* {i = 1,2,3) of class for some < a < 1 

with boundary (9F*. These are supposed to have a common boundary 

drl = drl = drl{=: i:,) (2.1) 



and fulfill the angle conditions from (1.3). As above, we introduce notation such that the 
outer conormals z/* at dTl fulfill 

l'i^l+l^i^! + l'i^! = on , 

and the normals A^* of F* are chosen such that 

^'Nl + j^N^ + j'N^ = on . (2.2) 



Note that we do not assume F* to be a stationary solution of (1.3), that is the mean 
curvature of F* can be arbitrary. 
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Let : f2* — > ]R"+^ be a local parametrization with C F* where is either an 

open subset of M" or -B^(O) = {x G M" | |a;| < Ir , > 0} in the case that we parametrize 
around a boundary point. For a e F^, we set F~^{a) — {xi{a), . . . ,a;„((T)) e R". Here 
and hereafter, for simplicity, we use the notation 

w{a) =w{xi,...,Xn) (cxeF*), 

i.e. we omit the parametrization. In particular, we set diw dxi{w o F). 
To parametrize a hypersurface close to F* , we define the mapping through 

¥ : Ft X (-£, s) X {-S, S) R"+' , (2.3) 
(a,w,r) 1-^ '^'(a,w,r) := a + wNl{a) +rrl{a) , 

where r* is a tangential vector field on F* with support in a neighbourhood of 9F* , which 
equals the conormal at c?F* . The index i has range 1, 2, 3. 
For i — 1, 2, 3 and functions 

: F: X [0, T] ^ (-£, e), : x [0, T] ^ {-S, S) 

we define the mappings = often omit the subscript {p^^iJi^) for shortness) 

through 

$^ F: X [0,T] ^ M"+^ , $^(a,t) := ^X^, /^'(pr'(^), ^)) ■ 

Herein pr* : F^ — >■ 9F* is defined such that pr* (cr) G dV\ is the point on dV\ with shortest 
distance on F* to a. We remark here that pr* is well-defined and smooth close to dT\. 
Note that we need this mapping just in a (small) neighbourhood of dV\. because it is 
used in the product /x'(pr'(cr), t) r*(cr), where the second term is zero outside a (small) 
neighbourhood of 9F^. For small £, 5 > and fixed t we set 

and finally we define new hypersurfaces through 

F^,,^,(t) := image(($Ot) ■ (2-4) 

We observe that for = and = the resulting surface is simply Fpz=o,;uiso(^) — 
for every t. 

Remark 2.1. We remark that for G C'^iV^) and jJ' G C^(E^) small enough in the 
C^{V\)- resp. C^{T:^)-norm the mapping ($*)* is a local C'^ -dijfeomorphism onto its image. 

In fact, omitting the time variable t and the index i for the moment, choosing a local 
parametrization and using the above abbreviations we calculate 

di^ = diq + dipN^ + pdiN^ + di{po pr) n + {po pr)din ■ 

A rather lengthy, but elementary calculation for gik = 9^$) gives 

gik = {g*)ik + Pik{P: ilJ' o pr), Vp, V(p o pr)) , 
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where Pik is a polynomial with Pik{0) = 0. With the help of the Leibniz formula for the 
determinant we can then derive 



g = det {{g, 



Ik Lk=l,...,n 



g* + P{p, (/u o pr), Vp, V(/i o pr)) 



where P is a polynomial with P{0) = 0. Since g^, > we conclude that for p and fi small 
enough in the C^-norms also g is positive. Together with the fact that {gik)i,k=i,...,n is 
positive semi-definite due to 



J2 ^iSikik = ^lidi'^^dkHk) 

l,k=l 



Lk=l 



1=1 



> for all ^ e 



(2.5) 



we conclude the property that {gik)i,k=i....,n is even positive definite. Hence we obtain a 
strict inequality in (2.5), whenever C, ^ arid we conclude that . . . are linearly 

independent, which means that the differential d^{a) has full rank. 

Finally with the help of the inverse function theorem we conclude that ($*)t is a local 
diffeomorphism and the image Ti{t) has metric tensor {gik)i,k=i,...,n- 

In the definition of we allow at the triple junction for a movement in normal and 
tangential direction, and hence there are enough degrees of freedom to formulate the 
condition, that the hypersurfaces Ti{t) meet in one triple junction S(t) at their boundary, 
through 

$i(a,t) = $2((T,t) = $='(a,t) foraeS, ,t>0. (2.6) 

We rewrite these equations in the following lemma, which was shown in Depner and 
Garcke El- 



Lemma 2.2. Equivalent to the equations (|2.6|) are the following conditions 

on , 



(i) + 7^P^ + 7^P^ 



[It) n 



[c^ fp — c'^p^^ on S^, . 



(2.7) 



for {i,j, k) = (1, 2, 3), (2, 3, 1) and (3, 1, 2), and where s* = sin^* and c' = cos 6^ 
With the notation ^ = {p^,p^,p^), p = {p^,p'^,p^) and the matrix 

^2 



r 







r.3 



v 







we can state the linear dependence from (ii) of (2.7) as 

p = Tp on . 



(2. 
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2.2 The nonlocal, nonlinear parabolic boundary value problem 



From now on, we always assume condition (2.6). We introduce the notation N^{a,t), 
V^{a,t) and H^{a,t) which are the normal, the normal velocity and the mean curvature 
of r*(t) := Tpi^^i{t) at the point $*((T, t). Then we write equation (1.3) over the fixed 
hypersurfaces Tl, Tl, and Tl as follows: 



V'{a,t) = (3'H'{a,t) 

{N^{a,t),N^{a,t)) =cose^ 
(pH^^,0),/i*(a,0)) = (p^,/i'o) 



on Tl,tE[0,T], 2 = 1,2,3 

on , t G [0, T] , 

on , t G [0, T] , 

on ri X , i = 1,2,3 , 



(2.9) 



where we assume that the initial surfaces Fq from (1.3) are given as 

F^ = image{a ^ ^'^(a, p^(a), /i^(pr*(a))) | a G F^}. 

Herein we assume Pq G C^+"(F*) with ||pq||(72+<:« < e for some £ > 0, pg G C'^^°'(J2^) given 
by = Tpo oil ^* cLud in addition the angle conditions from (1.3) for Fq shall be fulfilled. 
Furthermore, we assume that 



-f^/3^H^ + i^P^Hl + i^P^Hl = on 



(2.10) 



where Hq is the mean curvature of Fq. 



Note that equation (2.10) follows for smooth 
solutions from the first line in problem (2.9) at t = on S*, since for points on the 
triple junction we can write for the normal velocity = (c'(0),A^') with one curve 
c : [0,to] M"^^ on with c{t) G S()f:) and use equation (1.6) for Fq which follows from 



the angle conditions. 

Remark 2.3. The requirement that the C'^~^°'-norm of the initial values pg is small implies 
that the initial hypersurfaces Fq are 0"^^°" -close to the reference hypersurfaces T\, which 
are of class C^"*"". In order to make this compatible to condition (2.10), there are two 
possibilities. 

On the one hand we could start with initial hypersurfaces Fq, which fulfill (2.10) and 
then choose hypersurfaces T\ of class C^^'^ , which are close enough to Fq. This would 
imply that condition (2.10) is almost fulfilled for T\ in the sense that |X]^=i 7*/^*-^*| ^-^ 
small. 

On the other hand we could additionally require condition (2.10) for the reference 
hypersurfaces T\. In this way the above approach would always work in the sense that 
there are hypersurfaces Fq given by Pq, such that ||pQ||c2+ti is small and (2.10) holds. 

Due to the condition 9^ + 9"^ + 9^ = 2ti and the fact that the surfaces all meet at a 
triple junction at their boundary, which follows from (2.6), the third angle condition 

(iV2(a,t),iV3(a,t)) = cos ^2 , t G [0,T] , 



(2.11) 



is automatically fulfilled and we omit it from now on. The equations (2.9) give a second 
order system of partial differential equations for the functions (p^, p^, p^, p^, p^, p^). 



8 



More precisely, we can obtain the following representation for the equation. For the 
normal velocities we calculate 

V'{a,t)^{N'{a,t),dt¥{a,t)) 

= {N'{a, t), dt{a + p\a, t) Ni{a) + f^'{pAa),t) rtia)}) 

= (7V^(a,i),A^:(a)) + Wa,i),T:(a)) V(pr'(<^),^) ■ 

We remark that there is a function N'' such that 

N'{a, t) := N'{a, p\a, t), i,'{pT\a),t), Vp\a, i), W(pr'(^), t)) 

is the unit normal vector field of r*(t), where Vp* is the gradient of p* on the hypersurfaccs 
r* , which is denoted in a local chart by Vjp* = djp^ {j = 1, . . . ,n), and V/^* is the (n — 1)- 
dimensional gradient of on a surface E*. A formula for TV* can be given with the help of 
a^local chart through a normalized cross product of the tangential vectors 9/$'. Therefore 
A^* is a nonlocal operator, since in its formula wc find an expression /x*(pr*((T), t)r^((j) 
so that we do not only need p, ^t^and its derivatives at the point a but also the point 
pr((T) e dVl in order to calculate N\ 
Since 

{g% = {dj¥, dk^'), {h%k = {N\ djdk^'), 
the mean curvature if* = {g^y^ {h'') jk is represented as 

where V^p* is the Hessian of p* on hypersurfaccs F* defined in a local chart by 

= V,,V,,p* = dAP' - ^iiAP' = 1, . . . 

where F^^^^ are the Christoffel symbols for F* and we used the sum convention for the last 

2 

term. The expression V p.* denotes the Hessian of p.* on the (n — l)-dimensional surface 
E*. Note that the coefficients in front of the term V^Vfc p^ in are given by 

Thus the mean curvature flow equations can be reformulated as 

dtp' = a'{a, p\ p')H\a, p\ p') + ai(a, p\ (2.12) 

where H^{a, p\ p'^) := H'{a, p\ p\V p\V p\V'^ p\v'^ p"^) and 

/3* 



a\a,p\p'):=a\a,p\p\Vp\Vp') 
P\ P") ■= m<^, P\P\ Vp*, Vp^) 



{Ni{a),Ni{a,p\p\Vp\Vp')) 
{Tt{a),N\a,p\p\Vp\Vp')) 



{Ni{a),Ni(a,p\p\Vp\Vp')) 
Note that we omitted the mapping pr* in the functions p' for reasons of shortness. 
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Now we will write equation (2.12) as an evolution equation, which is nonlocal in space, 

oy u 

(2.13) 



solely for the mappings by using the linear dependence (2.8) on S^,. To this end, we 

and rewrite (2.12) into 



use (2.8) in the form /x* 



dtp' = r{p\ pIsJ + a]{p\ p\^Jdt{Tp o pr 

where (omitting the t-variable for the moment) 

r{p\p\^J{a)=a\a,p',{Tp\^y)H\a,p\{Tp\^y) 
a]{p\p\^J{a)=a]{a,p\{rp\^J) 

With the following notations on T,^, given by 

Tip, p|^J(a) = (^Xp\PlEjM).=l,2,3 

Pt(p, pIsJ(^) = diag [{a]ip\ p|sJM),=i,2 



for a e Ti 
for a G r: 

for cr G , 
for cr G 



(2.14) 



(2.15) 



we can write (2.13) as vector identity on through 

dtp = J^{p, pIsJ + V^{p, p\^JT{dtp) . 

Rearranging leads to 

{Id - V^{p, pIsJT) dtp = Tip, pIsJ on . 
Then, with the help of V{p, ) given by 

vip, p\^j := r {Id - i?t(p, Pi^jry' , 

it follows that 

TdtP = V{p, p\^jT{p, p\j,J on . 
In a neighbourhood of S^,, where pr* is defined, this leads to 

dtp\pT\a)) = {Tdtp{vr\a))y = ({P(p, p\^jnP. pIeJ} °PrT • 
Hence, the equation ( |2.12[ ) is rewritten as 

dtp' = r{p\p\^j + a'^{p\p\^:){{v{p,p\^jHp.pkJ]ov^Y on r:. 

The second term of the right hand side of this equation contains non-local terms including 
the highest order derivatives, that is, the second order derivatives. 
The angle conditions at the triple junction S,,, can be written as 

g\p):={U\p)M\p)) -cose'' = Q on S,, t>0, 

G\p) := {^^\p),^^%p)) - cos^^ = o on , t > o 

with the notation TV' (v) (a, t) := N\a, v\a, t), (r(i;opr(a, t))y , S/v\a, t), V(r(i;opr(a, t))Y). 
Note that due to a = pr*((T) for cr G the operators and Q'^ are local differential 
operators and depends only on p^ and p^ as well as only on p^ and p^. 
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Finally we have to take care of the equations (2.7), which are needed to make sure that 



that we are left with [2.7 given by 



the attachment condition (2.6) holds. Equation (2.7)(u) is already included implicitly, so 



g\p) ■= -i^p" + 7V + 7V = on , t > . 
Altogether this leads to the following nonlinear, nonlocal problem for i = 1,2, 3: 

dtp' = rip\ p\^j + a]ip\ p\^j {{Pip, pi^jTip, p\^j} opT^Y oil r:, t>0, 

g\p) =0 on , t>0, 

p%,o) = pi on r:. 

(2.16) 

2.3 The compatibility conditions 

For Pq we assume the compatibility conditions 

3 



g\po) = on and J] ^^^'^-^o. PoIeJ = on 



(2.17) 



i=l 



where /C* denotes the right side of the first line in (2.16). To state all the dependencies 



explicitly, we remark that by construction there is a function /C* such that 

Vip\ p\^Jia,t) = V{a,p\a,t),Vp\a,t),V'p\a,t), p\^^ (pr*(a), t), . . . 

V Pis, (pr^(a),t), Pis, ipr\a),t)) . (2.18) 

Note that we always set /Xq = Tpg on T,^ and therefore the geometric compatibility 



condition (2.10) is fulfilled since we require (2.17) for pg. This is stated in the following 
lemma. 



Lemma 2.4. The compatibility conditions (2.17) for pg imply the geometric compatibility 
condition (2.10). 

Proof. Using the abbreviations /Cq = 10{pQ, PolsJ and £q = (T/Co)\ where /Co = 



(/Co)i=i2 3' from the second compatibility condition in ( 2.17[ ) with arguments sim- 

ilar as in the proof of Lemma 2.2 (see [8]) that 



lC\,Nl + C,Tl = lClNi+iy,rl on . 
Now we show on S^, the following identity 

{{1C^NI + C,T^,NI) = I3'HI on . 



(2.19) 



To see this, we write in the following an index on every term to indicate evaluation at 
Po to get 

= ai^w, + ai, {r{id - v^^ory'PoY 

/Co = J-Q + ViflTild - Pt.oT)-^ J-Q , 
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respectively. With the definition of a* and aj this leads to 

/c^(iv:, iv^) = p^Hi, - (r:, iv^) {T{id - v^^oTr'ToY . 

In order to obtain ( |2.19[ ) it is therefore enough to show that 
which is, without loss of generality, equivalent to 

r(/d-Dt,or)-^j-o = r/Co. 

To obtain the last equality we observe that 

{Id - I^t.oT)-^ J-Q - /Co = {Id - I?t,or)-^ J-Q - J-Q - 2^t,or(/rf - l^t-oT)-' J-Q 

= {Id - v^^,r){id - Pt,or)-^j-o - J-Q 

= 0, 



so that finally (2.19) is verified. 

Since the term in brackets on the left side of (2.19) is independent of i, we can multiply 
by 7*, sum over i = 1,2,3 and use (2.2 ) resulting from the angle conditions for Fq to derive 
finally equation (|2.10D, that is Y.]=i YP'H'q = on S,. □ 



3 Linearization 



In this section we will derive the linearization of the nonlinear nonlocal problem (2.16) 
around p = 0, that is around the fixed reference hypersurfaces F^. This will be done by 



considering the geometric problem (2.9 ) and linearize this around (p, p) = 0. For this part 
we can use the work of Depner and Garcke [8], where the authors considered stationary 
reference hypersurfaces, and comment on the differences. To explain our notation we give 
the calculations for the normal velocity and just refer for the linearization of the mean 
curvature and the angle conditions to j8]|. In each term in (2.9), we write eu^ and £0* 
instead of and /x* for i = 1,2,3, differentiate with respect to e, and set e = in the 



resulting equations. Here, we have to assume the triple junction condition (2.6) for ^i, 
which is nothing else than assuming it for . In this way, we will get linear partial 

differential equations, where we then express terms of 0* as nonlocal terms in u with the 



help of (2.8) for u and <p. 



Linearization of the normal velocity: For the linearization of the normal velocity 
we obtain 



d 

de 



£ = 



{Nl{a),Nl{a)) dtu\a,t) + (iV:(a), r:(a)) dt<P\pT\a),t) 



e=0 



dtu\a,t). 
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Linearization of the mean curvature: For the linearization of the mean curvature 
if*, we use the following result, see Depner, Garcke [8j and Depner [7j, where [7j contains 
the detailed calculation: 



d 

de 



e=Q 

Ti 1 2 



de 



e=0. 



) 



where Ap^ is the Laplace-Beltrami operator on F*, 11* denotes the second fundamental 
form of F* and |n*p is the squared norm of 11* and hence given as the sum of the 
squared principal curvatures. Furthermore Vr* is the surface gradient on F* and [■]'^ 
is the tangential part of a vector. Note that the last term would vanish for reference 
hypersurfaces with constant mean curvature. For the last term we compute 



de '""''"f 





T 


e=0. 





d 

de 



[a + eu\a, t) Nl{a) + £0*(pr*(a), t) rl^a)) 



e=OJ 



: [n*(a,t)iV:(a) + 0*(pr*(a),t)r:(o-)]^ 
:0*(pr*(a),t)r:(a). 



so that we get 

d 
de 



if*o$*^._^^.(a,t) 



£=0 



Ar.M*(a,t) + \ni\\a,t)u\a,t) + lyriH\a),Tl{a)) 0^(pr*(a),t) . 



Linearization of the angle conditions: The linearization of the angle condition 
(A^*, A^-') = cosO'' is the technically most challenging part and we use the following result 
of Depner and Garcke [8] : 



d 

de 



/ AT* o $* ■ , ■ A^-^' o . ■) 



£=0 



= d^.n' + n:(z/, z/)0* - - Uiiui, vlW 

on for t > and for = (1,2) and (2,3). Note that in [8j there was a second 

equivalent formulation of the above formula, which is not possible here, since the reference 



hypersurfaces are not stationary. Nevertheless with the help of (2.8) we can get rid of 
by expressing it with the help of u. 



Altogether, we get for the linearization of (2.9) the following linear system of partial 
differential equations for («*, 0*) and z = 1, 2, 3. 



' 5tu* = /3* (Ar>M* + in* pM*) + /3*(Vr.if*, r*) (0* o pr* 



'^^u" + 7^u^ + 'X'u^ = 



{u\ 0* 



li=0 



(Po./^o) 



on F* X [0,T] , 
on X [0,T] , 
on X [0,T] , 
on X [0,T] , 
on F*. 



(3.1) 
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Note that 0* o pr* can be rewritten as {l~{u o pr*))* due to equation (2.8), which also has 
to hold for and u. Now we are able to rewrite the nonlinear, nonlocal problem (2.16) as 
a perturbation of a linearized problem. Let the operator and the function Q be given 
by 

A' = /3*{Ar. + |n:p/}, C(a) = /3^(Vrii^:(a), r:(a)) . 

We also introduce an operator corresponding to the linearized boundary conditions given 
by 



K 



'cm — c'u* 



where := H^z/^, i/*) denotes the normal curvature of V\ in direction of v\. 

With this notation we can rewrite the nonlinear nonlocal problem (2.16) into the 
following one, where i = 1, 2, 3: 



' dtu' = A'u' + C{r{u o pr^))^ + f{u\ u\^J on x [0,T] 
3 



u 



\t=o ~ Po 

Herein, f and b* are defined through 



on X [0,T], 

on ri. 



(3.2) 



r{v\ v\^J := r{v\ v\^J-\ AV + C{T{v o pr^)) 



bHv) :-- 



+ a\{v\ v\^J{{r{v, v\^jT{v, v\^J} o pry 
3 



(3.3) 
(3.4) 



Note that the first boundary condition on the triple junction S,,, in problem (2.16) is 
already linear and therefore b^{v) = 0. But we will nevertheless use to avoid some case 
by case analysis. 



4 Analysis of the linearized problem 



In this section we consider the linear nonhomogeneous problem corresponding to (3.2). 
We will give a local existence result for the case with initial data zero and then outline 
the necessary steps for the arbitrary case. First we introduce for an arbitrary smooth 
Riemannian manifold {T,g) some notation. For an integer k and smooth functions u : 
r — > M, we denote by V^u the k-th covariant derivative of u and by \'V'^u\ the norm of 
V'^u defined in a local chart by, see e.g. [1], 
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Note that VjU = diU and Vj-^^i^u = Vi-^Vi^u = di-^di^u — V^^^dmU. For T > and 
< a < 1, set Qt = r X [0, T] and 



sup |M(cr, t)|, 



sup 



[^((j, t) — u{a,t)\ 



sup 

(o-,i), (a,t)eQT,t=^t 



\t-t 



where dg denotes the distance on F induced by the metric g. Then, we define the norms 



\u\ 



and ||M||p2+.,i+f (Q^) 



as 



\u\ 



\u\ 



t 1 



c 



2 + a,l + 



+ llVulL + llV^ul 



+ 



Set Xt = C2+"'i+t(g^) X C2+-'i+f (g2,) x C2+«'i+f (g3,), where = F: x [0,T]. 
Then we have the following theorem about existence of solutions to the linearized, non- 
homogeneous problem with initial data zero. 

Theorem 4.1. Let a G (0,1). Then there exists a 5q > Q such that for every /* G 
C"'^(<5U ^ Ci+"'H^(S, X [0,5o]), i = 1,2,3, with = and which fulfill the 

compatibility condition 



il'f' + l'P + l'f') 



0, b' 



t=o 



on 



2,3, 



the problem 



dtu' = A'u' + C(T(n o pr^))^ + f on Fj, x [0, T], 







on X [0,r], 



on F! 



(4.1) 



for i = 1,2,3 has a unique solution {u^,u^,u^) G Xs^. Moreover, there exists a C > 0, 
which is independent of Sq, such that 



i=l 



i=l 



1 I i+Q r . 



First, we will consider problem (4.1) without the nonlocal term (T('U o pr*))* and at 
the end we will include it with the help of a perturbation argument. 

In order to apply the C"-regularity theory of Solonnikov |28] we need to show that the 
boundary value problem (4.1) fulfills the Lopantinskii-Shapiro compatibility conditions, 
see Chapter I of [28j, where the conditions are stated. To this end we have to rewrite 
problem (4.1) with the help of local coordinates and a partition of unity as a problem 
in Euclidean space. We will do this locally around the triple junction with specifically 
chosen local coordinates, since the compatibility conditions have to be checked just there. 
Locally around a point cr G we choose for each of the surfaces F*, i = 1,2,3, local 
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coordinates (xi, . . . , x„) such that {xi, . . . ,Xn-i) parametrize S^, and such that the metric 
tensors fulfill 



{9')^^ = l, (^7%„ = for j = l,...,n-l 



(4.2) 



This is possible by choosing the n'th coordinate as the distance from the (n— l)-dimensional 
surface S^,. 

Denoting the representation of the u\ j = 1, 2, 3, in local coordinates as u\ j = 1, 2, 3, 
the principal parts of the boundary operators in (4.1) can be written as 



BqUj = < dnU^ - dnU^ , 
i=l I dnU^ - dnU^ , 



i = l, 
i = 2, 
i = 3. 



The principal part of the parabolic differential operator takes the form 

£o(9„V) = (/;0,,=,,,3 

with 



0, t ^ j, 

n 

dt-Y,P'g''''d,di, i=j. 



kl=l 



For ^ G and p G C with positive real part we now define 

3 



L := det£o(p,ie) = \{{p+Y. P'9"''^k^i 

i=l k.l=l 



and 



Lemma 4.2. The operators {Co,Bq) fulfill the Lopantinskii- Shapiro conditions. 



Proof. For the coefficients of Co we calculate 

0, 



7 U 

'■0 



« 7^ J, 



6 n 



j=i 



k,l=l 



We now set ^ = ^' + re„ with = 0, r G M and e„ = (0, . . . ,0, 1). Let t\p,^'), i = 1,2,3, 
be those roots of L{p,i{^' + re„)), which have positive imaginary part. The fact that 
there are exact three roots with positive imaginary part follows from the fact that the 
system in the first line of (4.1) is parabolic. Now we define 

n-l 



kJ=l 
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where p is assumed to have a positive real part. We choose polar coordinates 

The fact that p has positive real part and the fact that i^g^'^''^ is positive definite 
imply that 0' G (— 7r/2, 7r/2) . Hence we compute 

r'(p,0 = V^y^e^^(*'+-). (4.3) 

The Lopantinskii-Shapiro conditions now require that the rows of the matrix BqCq are 
linearly independent for all p e C with Rep > modulo the polynomial 

3 

M-'{p,e,r) = l[{T-T\p,0} ■ 

4 = 1 

This can only be true if 

3 

J2^'Bi^{r)=0 mod r - r^(p, ^ , J = 1, 2, 3 

i=l 

has a nontrivial solution {co^ , co"^ , co^) , where Bq {j) = BQ^ren)- Hence we need to decide 
whether the set of equations 

Y,^^Bi,\T\p,a) = o (4.4) 
1=1 

has a nontrivial solution. Using the definition of the we finally need to decide whether 
the determinant of the matrix 



7^ 
7^ -T^ 

is singular or not. Here we abbreviated r* = T^{p,^')- The determinant is given as 

7W^ + 7W3 + 7W^ (4.5) 

In polar coordinates the angle of rV-' is given as {(j)^ +((>') /2+7i. Since 0*, 0^ G (— 7r/2, n/2j 
we obtain that rV-' has negative real part. Hence 7^r^r^ + 7^r"^r^ + 7'^r^r^ is the sum 
of three summands which all have negative real part. Hence the determinant is non-zero 
and we have shown that the Lopantinskii-Shapiro conditions hold. □ 

Remark 4.3. In Latushkin, Priiss and Schnauhelt f2^ the Lopatinskii- Shapiro condition 
is formulated as a condition for a system of ordinary differential equations. In our notation 
this reads as follows. Let a G and (xi, . . . ,x„) be local coordinates (in a region Q) as 
in (4.2) and set ^o(V) = diag((— X]L=i l^^ 9^'^^^ ^kdi) 1=1,2,3) ■ Then the formulation in 
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requires that for given C, € MJ^ with ^ -L n and A G {-2 G C | Re{z) > 0} with {X,^) 7^ (0, 0) 
the function (p = is the only bounded solution in Co(M+; C^) of the ODE-system 



\^{y) + + n{x)dy)^{y) = ^ , y>0, 



(4.6) 
(4.7) 



The equivalence of the formulation in to the algebraic formulation in Solonnikov 
can be found in Eidelman and Zhitarashu [TU. Chap. 1.2]. 



By choosing for simplicity as above ^ = (^',0) and n{x) = Cn the equations (4.6) 
and (4.7) reduces in our case to 



(4.8) 

= 0, y = 0, (4.9) 

(/)' = (¥^7 = (y^'y, y = 0- (4.10) 

These equations can be treated with an energy method to show that a solution must be 
zero. To this end we test line (4.8) with j^Ip^/P^ and sum over j = 1,2,3 to get 

= 5^(A + /3^irP)^ / y\^dy-Y,l' i^r^'dy 
,=1 Jo ^.=1 Jo 

dy+Y.^' / w)'?dy-Y,i'iv')'mK^) 

Jo 



3 3 r'^ 



3 



^ poo '-' 

dy + J2i' / |(y.7l'rf2/-(^'y(0)5^7V(0). 



/n i/ie /as^ /me we used the boundary condition (4.10). Finally with (4.9) we see that the 
last term vanishes and that therefore {ip^ , ip"^ , (f^) = 0. 



Proof of Theorem 4-i First we construct a weak solution of problem (4.1) without the 
nonlocal term. In order to apply an energy method we modify the equations into 



,dtu' = 7*(Ar.M' + in^pM*) + -^f on x [0, T] 



1 
3 



r 



u 



t=0 







on X [0,T] 



on n . 



(4.11) 



In this way we are able to choose the weak solution u = {u^, u^, u^) and the test functions 
^ = (^^5^^;^^) ill the same space. Now we introduce the function spaces 

£ := L\Tl) X L\Tl) x ^^(r^), := ^^(Sri) x L\dTl) x ^^(Sr^) (= (^^(S,))'), 
:= H\Tl) X /Ji(r,2) ^ i/i(r3), ^ := {t^ g | + 7V + 7V = a.e. on S J. 

Also, we introduce the time-dependent bilinear form 
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for u{- , t) ,$,{■, t) G S. The weak formulation then reads as follows. Find u G L'^{0,T]S) 
with dtu G L2(0,T; (H^)"^) such that 

{dtu, Oduai + B[u, ^; t] = (/, + 6(1; t) for all ^ G £: and a.e. in t, (4.12) 
where (T-L^) ^ is the dual space to "H^ and 



{dtU,Odual = if,^)c = 5^^(f ,e)L^{n) (4.13) 

i=l ^ i=l ^ 

are scaled versions of the corresponding duality pairing and inner product. The time- 
dependent linear form b is given through 



and consists of terms which appear formally due to the rewriting of J^, Y^ui^'^C d'H"' ^ 
to make use of Yl^=i 7*^* ~ 0- That this weak formulation for smooth solutions is equiv- 
alent to the strong formulation, can be checked by a straightforward computation using 
integration by parts and the restriction ^ G 

We want to apply the Galerkin method and therefore assume that Wk = Wk{o') for 
/c = 1, 2, . . . are smooth functions such that {lOfcjfcLi is an orthonormal basis in C Indeed, 
we can take such {lUfc}^^ considering eigenfunctions of the eigenvalue problem 

-yAp. = A '^y on r: , 2 = 1, 2, 3 , 

P 

7^w^ + 7^w^ + 7^w^ = on S^, , 

{Vviw\ vl) = {Vr^y, vl) = (V^3«7^ i^!) on . 

This follows similar as in Gilbarg and Trudinger [iHj by considering the quadratic form 

3 



Qiu,u) = J2Y [ {Vr^u\VT^.u')dK 
i=i -^n 



on S and the norm of C as in (4.13). In addition the eigenfunctions are orthogonal 
with respect to the quadratic form Q. We remark, that since the boundary conditions 
fulfill the Lopantinskii-Shapiro conditions, one can also derive regularity results for the 
eigenfunctions {lo^}^^. 

Now fix a positive integer m G N and look for : [0, T] — t- £^ of the form 

m 

u^{t) = Y,<t{t)wk. (4.14) 

fe=i 

Here the coefficients (i^(t) for k = 1,2, ... ,m have to be chosen such that 

dtiQ) = 0, (4.15) 

{dtUm, Wk)c + B[Um, Wk] t] = (/, Wk)c + b{'Wk] t) , (4.16) 
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where k = 1, . . . ,m and the second hne has to be understood pointwise in t. Note that 
due to Wk G £ a function Um of the form ( |4.14[ ) satisfies 

j'ulia, t) + j^Ua, t) + i^vt{a, t) = for a G . 

With the help of theory for hnear systems of ordinary differential equations we find 
(li^, . . . , d™) as a unique solution of 



(4)' W + ^[^'' ^]4(^) = ' ^)^'^k)c + h{wk- 1) 



1=1 



with the initial data (4.15), so that Um of the form (4.14) satisfies (4.15) and (4.16) for 
each m G N. 

Since the trace operator is compact one can use a contradiction argument similar as 
in the proof of the Ehrling Lemma in order to derive the inequality 



\u\ 



< e\Wu\\l + CJu 



|2 



for each e > and a constant > 0. Using this inequality one can argue similar as in 
the proof of Evans |T2l Sect. 7.1.2, Th. 2] and obtain the energy estimate 



sup + \\Um\\L^{0,T;S) + \\dtUm\\ L^(o T-m^)-^) - ^ {\\ f\\ L^{0,T:C)) + II ^11 L2(0,T;Z;i,) ) 

0<t<T 

(4.17) 

for m G N and a constant C > 0. Using this we can prove the existence and uniqueness 
of a weak solution with standard arguments, which can be found for example in Evans 
[m p.356-358]. 



Let us derive Schauder estimates for solutions of problem (4.11). Here we consider the 



Holder estimate only near the triple junction and just remark that away from the triple 
junction the result follows in a standard way after localization. 

Let us introduce some notation. Locally around a point cr G we choose parametriza- 
tions which flatten the boundary in the following way. We pick a sequence < ri < r2 < 
rg < r4 and with Qi := Br,{y) n {x G M" | x„ > 0} for / = 1,2,3,4, where y e W 
is such that yn = 0, we let : (^4 — )■ L*, i = 1,2,3, be local parametrizations with 



F%y) 



a and F' 



{Xn=0} 



C S*. Additionally for a given to > we choose a sequence 



< 5i < ^2 < ^3 < ^4 and set A, := (to - Si,to + 6i) f] {t e R\t > 0} for / = 1, 2, 3, 4. 



With the help of a cut-off function we will formulate problem (4.11) for the repre- 
sentations = o in X A4 in Euclidean space. To preserve the structure of the 
problem and to keep the notation simple, we will identify the notation of the function 
with its representation in local coordinates. In the next steps the sets Qi x A^ will be 
successively reduced to achieve finally the stated Holder estimate in Qi x Ai. We will 
need the following notation for parts of the boundary of Qi: 

Ci:=dQir\{xeW\xn>Q} and Si := dQi\Ci . 

Now let ?7 be a cut-off function satisfying 

r/ G Co°"(Q4 X A4), 0<r7<l, = 1 on Q3 x A3 . 
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We remark that due to the fact that (^4 is not open, the values r]{x,t) for x & S4 do not 
necessarily vanish. The same holds true for ri{x, 0), if A4 is not open. 



Now set = rju^, where {u ,u ,u ) is a weak solution of (4.11) and note that we do 



not distinguish between the functions -u* and its representations. Then we have in a weak 
sense 



Since (m^, u^, u^) is a weak solution of ( |4.11[ ), we deduce that {v^ , f ^, v^) is a weak solution 
of 



X 



dtv' = -f%A^^v' + |n:| V) + f{x, t), 



e Q4 X A4, 

0, {x,t) G 5*4 X A4, 

{Vriv\ vl) + nlij vf - {Vviv\ vl) - kI{T vf = b\x,t), {x,t) e x A4, 

(V^2^;^ u!) + kI{T vy - (V^3^;^ u!) - kI{T vf = h\x, t), (a;, t) G ^4 x A4, 
v\x,t) = 0, {x,t) G C4 X A4, 

^ 0) = 0, X G Qi, 



(4.18) 



where i = 1, 2, 3 and 



Note that 



/iQ3xA3 = f eC"^'t(Q3xA3), f 



Q4XA4 



G L\Q^ X A4) 



l + Q 



&I53XA3 = e C^+"'^(S3 X A3), b^s.xA, G ^'(54 X A4) 

Let and 6^ be smooth approximations of /* and 6* satisfying 
and on Q2 x A2 C (^3 x A3 we require 



(4.19) 



II/, 



< 



II&: 



C"'3r{Q3xA3) 



1 + a, 



^(S'sXAg) 



C"'7(Q3XA3)' 



C 



l + a. 



'iSaxAaY 



Replace /* and 6* by and 6^ in (4.18), and call this problem (4.18 )„. Since we checked 



the Lopatinskii-Shapiro conditions on the triple junction in Lemma 4.2 we can apply 
results from Solonnikov [28], Theorem 4.9] to get a unique solution G C'^~^°''^~^'^ {Q4 x A4) 
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of problem (4.18 )„. Using furthermore the local estimate from [28, Theorem 4.11], we 
obtain for Qi x Ai C (^2 x ^2 C Q3 x A3 from above that 



E 

i=l 



\Vj^\\„2 + aA + 



(■ 3 ^ 3 ^ 

ixAi) - 5Z "•^"llc^'f (Q2XA2) + 5Z Il^"llc^+">T^(52XA2) 

1=1 i=2 



+ ^ IKI|l2(Q2xA2) 
i=l 
3 



< 



+E 



C°-7(Q3XA3) ' ^ II" llci+'^-^(S3xA3) 
1=1 i=2 



(4.20) 



+ C'2 ^ |K||l2(Q4XA4)- 

i=l 



By means of (4.19) and the energy estimate (4.17) for the approximated problem (4.18), 

we see 



E 

2=1 



\Vn\\L2(A^ 



,H1(Q4)) ^ L2(Q4XA4) n ||L2(54xA4) 

1=1 



< 



5 E 



i=2 
3 



L2(Q4xA4) + ^ II" ||L2(S4xA4) 

=1 1=2 
^3 3 

^ II/1l2(Q4xA4) 



i=l 



+ $^ll« II L2(54XA4) 
1=2 



(4.21) 



+ X] ||MlL2(A4,i?l(Q4)) 
■t=l 

3 3 

< C'< ^ ||/1l2(Q4xA4) + X] l|0 ||L2(54XA4) 



1=1 



i=2 



In the last inequality we used the energy estimate (4.17). From the last bound we deduce 
the existence of a subsequence {v^^} C {f^} and of E L"^ {K^, {Q 4)) such that 

^n, weakly, 



and (f^,f^,f^) is a weak solution of (4.18). By uniqueness of the weak solution of (4.18), 

= in Q4 X A4. 



Let us rewrite as f^. By (4.20) and (4.21), we obtain 



E 



1=1 



|^2 + a.l+t(Q 



.xAo < ci{y: 
^ i=i 

^ i=l 



+E 



C"^^(Q3XA3) ' ^ II" Ilci+"'"^"(53XA3) 
1=2 



3 3 

L2(Q4XA4) + II" ||L2(S4XA4) 

1=2 
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< 



IC"^7(Q4XA4) 



4 = 2 



^(54XA4) 



Then, by the theorem of Arzela-Ascoh, there exist {v} } C {v}} and -0* G C^'^iQi x Ai) 
such that 

vi^^v' in C^^\Q^xA^). 
Here -0* is in C^'^°''^~^^ {Qi x Ai) because of, for example, 

\V,Vkv\x) - V,Vkv\y)\ = hm |V,Vfct;i^(x) - V,VkvUy)\ < C\x - 

It follows from uniqueness of a limit and tJ* = f * in (^4 x A4 that 

■0* = tJ* = t;* in Qi X Ai. 
Since = in Qi x Ai, is in (7^"'""'-'^+? (Qi x Ai) and satisfies 



E 



ixAO ^ J] llf llc-?(Q,xA4) + Il^'llc^+-4^(54XA4) f • 
i=l i=2 ^ 



Hence we are led to the stated Holder estimate locally around the triple junction S^,. By 
a covering argument we can enlarge the estimate to a neighbourhood of S^, and then by 
an easier argument, that we omit here, we can give it for all hypersurfaces T\ as claimed. 

Finally, by a perturbation argument as in Baconneau and Lunardi [21 Thm. 2.3], we 
derive the existence of a unique solution and the Schauder estimate for the linearized 
system with nonlocal term. We omit the details since this part is even easier than in [2] 
due to the fact that the nonlocal terms {T{u o pr*))* do not contain derivatives of u. 

□ 



Altogether we proved Theorem 4.1 



Remark 4.4. For the case of arbitrary initial date u^\^_q = Pq, we have the following 
existence result. Let a G (0,1). Then there exists 6q > such that for every /* G 
C"'^ (QU; G (S, X [0, 5o]) with = and G C2+"(r:) with the compatibility 

condition 

3 3 



i=l 



the problem 



4 = A'u' + C(T(w o pr^))^ + f on V\ X [0, T], 



u 



t=o 



on X [0,T], 



on rl 



(4.22) 



for i = 1,2,3 has a unique solution {u^,u'^,u^) G Xs^. Moreover, there exists C > 0, 
which is independent of 5q, such that 



i=i " i=i 



+ IIPollc2+-(Qjj}- 



For the proof consider the difference f* := — pg and apply Theorem 4-1 



to v\ 
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5 Local existence 



With the help of the previous results we are now in a position to solve the nonlinear 
nonlocal problem (3.2) locally in time. We will apply a method similar to Lunardi [23| 
Th. 8.5.4] resp. Baconneau and Lunardi [2]. But since we do not linearize around the 
initial state and since our problem is geometrically more involved, we state some of the 
arguments in detail. Note that for T > and < a < 1 we use the Holder spaces 



3^ 

T) ' 



where Qlp = F* x [0, T]. Roughly we show in the following theorem that if the initial state 
satisfies the compatibility conditions and lies C^^"-close to the reference state, there is a 
unique solution (m^, m^, m^) G of (3.2) where 5 > is chosen sufficiently small. 

Theorem 5.1. Assume that G C^^"(r*), i = 1,2,3, fulfill the compatibility conditions 
(2.17). Then there exist constants Rq > and Eq > such that for each R > Rq there is 
6 > satisfying that if Y2^=i \\Po\\c^+°'iri) ^ ^O; the nonlinear nonlocal problem (3.2) has 
a unique solution u = in the ball Br{pq) C X^. 

Proof. Let r > be a constant such that for G C^(r*) with Yl^i=i Il'^'llc2(ri) < f the 
following assumptions hold: 



(Al) v\j. ) and a\{v\ v\j. ) (see (2.13)) are well-defined as well as V{v, ) (see 

(A2) Any first order derivatives of J-"* with respect to v\ ^j'^^i ^ i ''^Is.' VjVfcf* 

and VjVfc vj^ are locally Lipschitz continuous with respect to those. Also, any 
first order derivatives of a| with respect to f*, i^l^,; ^j"^* ^^"^ '^Ie locally 
Lipschitz continuous with respect to those. 

(A3) Any second order derivatives of b* with respect to v\^^ and Vj 'v\j.^ are locally 
Lipschitz continuous with respect to those. 



We remark that these properties are realized for sufficiently small r since with the 

, ... 2 

notations z\ = {v\ Vf*, v\^^ , V I'l^J a^id Z2 = {v\ Vf*, V^f\ v\j,^ , V v\^^ , V v\j.J the 
quantities {g'')jk, det (((7*)^^), A^*, and {h^)jk are represented as 



{9'),k = {gl),k + P;fc(4) , 9' = det {{g^),k) = det {{gl),k) + P^^l) , .^ 
= Nl R\z\) + Q\z\) , h], = (HD^k R\z\) + 5i,(4) , ^ • ^ 

where P*^ and P* are polynomial functions with P*fc(0) = and P*(0) = 0, and P\ and 



2.1 



we 



S*]^ are rational functions with P*(0) = 1, Q*(0) = and S^^j^ifi) = 0. From Remark 
know that (7* 7^ for f * small enough in the C^-norm and that . . . , are linearly 

independent, in particular 191$* x . . . x 7^ 0, and therefore also N'^ is well-defined. 

Now fix P > and define the set 

Vr = {{v\v\v^) G Xs I v\a,Q) = pl IK'* -Pollc7^+-i+f (Qj) < R}- (5-2) 
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For V G we deduce from a standard estimate for parabolic Holder spaces, see e.g. 
Lunardi [23, Lem. 5.1.1], that for all t G [0,T] we have 

3 3 3 



i=l 



i=l 



i=l 



< [5^ +C5— +5) \\v'-pl\\^2+. 



+ £0 



(5.3) 



1=1 



< C5^R + eo, 

where the positive constant C depends only on a and max{l, 6^^^}. This shows that for 
sufficiently small 6 and Eq the operators J-"', and b*, evaluated at functions of the form 
v^{-,t), satisfy (Al)-(A3) for all t G [0,6]. In particular we remark for later use that for 
the right hand side /C* of the first line in (2.16), which is a combination of terms of the 
form J-"* and a|, we can conclude an analogue statement as in (Al)-(A2). This means that 
for D, to G the operator /C* is well-defined and it holds 



1=1 



(5.4) 



where is any first order derivative in {d^i , , 8^2 ,,»,(9^,i , (9^ ..i ,3^2 . |. Note 
that L depends only on the chosen r > from the beginning of the proof. In particular 
the same estimate holds true for v = and w = 0, i.e. 



\D,)c\pi poIeJ - D.}c\o)\u <lY IIPoI 



c2+«(ri) 



(5.5) 



1=1 



Due to the Lipschitz-continuity we also have that D^IO is bounded as a mapping from 
C C^+"'^"^^((55) into C°''^{Q\), which will be used later to estimate 



nonhomogeneous problem for i = 1, 2, 3: 



Fix V = {v ,v ,v'^) G Vji and let u = {u ,u ,u 



[D^IC\v\v\^J]^^.^ <C{R). (5.6) 

= A(v) be the solution of the linear. 



' dtu' = A'u' + C{T{uopT')y + f{v\v\^J on rtx[0,5] 
3 

Y^''u^ = b\v) 



u 



t=0 



on T.^ X [0, 5], 

on ri. 



(5.7) 



Due to the compatibility condition (2.17) for Pq, we see that f and satisfy the necessary 
compatibility conditions to apply Remark 4.4, that is 

33 3 
Y^'f('^'^''\^X=o = -Y^'{^'Po + (i^Poy) and b^(t^)L=o = E^Vo on . 

1=1 i=l j=l 
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Therefore we get a unique solution -u G of (5.7) for given v G for a possibly smaller 
5 > 0, but not depending on the choice of v G Vj^. 

If we are now able to find a fixed point of A, then this is a local solution to the nonlinear 
problem (3.2). Thus we will prove that A maps Vji into itself and is a contraction for 
suitable 6, Eq and R. 

For V, it> G Pij, we see that u = A{v) — A{w) is the solution of 

' dtu' = A'u' + C{T{u o pr^))^ + f{v\ v\^J - f{w\ w\^J on Ti x [0, 5], 



= b\v) - b\w) 

i=i 

I n*(.,0) = 



on X [0,(5], (5.8) 



on ri 



for i = 1, 2, 3. Then, by means of Theorem 4.1, we have the estimate 



i=l 



3 



<C^[\\f{v\v\^J~f{w\w\^:, 



+ Wiv) - b'iw) 



l + a 



{S,x[0,5]) 



Now we claim that there are constants C {R) and L such that 



3 



i=l 



x[0,S]) 



(5.9) 



W „2 + a,l + 



where C{R) is independent of 5 and L is as in (5.4). To show the estimate for f , we use 
the notation = -A^v + CO'i'^ ° pr*))* for the linearization including the nonlocal 

terms to get, compare (3.3), 

Note that herein = 9/C*(0)i> is the linearization around the reference hypersurfaces 

represented through p = and that /C* is a nonlinear nonlocal operator depending on v\ 
Vf*, V^f*, v\j.^, V v\j.^ and V v\j.^, compare (2.18). 

The difference in f can be written locally with the help of a suitable parametrization 
as follows 



d 

w\ vis. , w\j.J) ds - Al^^iv - w) 



■■Q\v\w\v\^^,w\^J{v'-w^) + Y,Q 



k=l 
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3 n-1 



j=l k=l 

n 



k,l=l 
3 n-1 



^2 ' 7 1 



j=l k,l=l 

+ {dr{pip,\^j-dV{o)) {v~w), 

where with = {Po, PoIeJ following notation 

w\ v\^_^ , w\^J = {sv' + (1 - s)w\ s v\^^ + (1 - s) w\^J , 

&{v\w\ v\^^ , w\^J = [ {^^^lC{is{v\w\ v\^^ , w\^J) - ^,,^)C'\^o)) ds, 

Jo 

^i{v\w\v\^^,w\^J = [ {^^,y^lC'{^,{v\w\v\^^,w\^J) -ds/,^.IC'{^o)) ds, 

Jo 

Herein, by a slight abuse of notation, we identify the /C*-terms with its localized versions. 

Now we observe for 6 G {6*, 9*, 9^, 9^'"', 9^ 9^'^;} that Q\^^q = 0, and therefore we 
derive 

||9||oo<5^(9)?<C(i?)<5t. 



Additionally (5.5) gives 

3 3 

1=1 1=1 

so that we arrive at 

3 

^^ ^IsJ - f ^"IeJIIoo < (C(i?)5t + Leo) 11^' ~ ^ilc2+"-i+^(Qj) • 

i=l 

Moreover it follows from D{v' - w%^q = 0, where D e {V°, Vfc, V|;, V°, V^, J (of 
course for surface gradients V we restrict the function — to the triple junction E,,,), 
that 



i=l 
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Set [ • ](^c,,f = {■)x ~^ {')t ^iid Ist 9D(f * — w*) be corresponding to each other as in the 
formula for the difference in f . Then we obtain 



< C{R)62 ^ lit;* - W ||^2+a,l + 

i=l 



Additionally it follows from (5.5) and (5.6) that 



Thus we are led to 

3 



i=l 



\V — W 



i=l 



By using (A3) we can give analogously an estimate for the differences in b* and therefore 
we arrive at the inequality (5.9). Consequently, we obtain that A is a 1/2-contraction 
provided 6 and Eq are small enough. 

To see that A maps Vji into itself, we have for v G and u = A{v) 



1 1=1 

R 



< Yl {W^i'^y - ^(Po)1c2+<^.i+f (Q.) + l|A(Po)' - Pollc2+".i+f (Q.)) 



< 2 + Eii^(/^or 

i=l 



For the second inequality, we used the fact that A is a 1/2-contraction provided 6 and Eq 
are small enough. The function w = A(pq) — is the solution of 



f dtw^ = ^^ip^ + JC^ipl polsj on r: X [0, S] 

3 

t w*(-,0)=0 



on X [0,(5], 

on ri. 



(5.10) 



4.1 



are 



Due to the assumptions (2.17) on Po the compatibility conditions from Theorem 
fulfilled and we can apply it to get the existence of a C > independent of 5 > 0, such 
that the solution of (5.10) satisfies 

3 3 



\W ||^2 + a,l+t(Q«) 



i=l i=l 

We estimate the right side of the above inequality by C = C'{eo) and we arrive at 

R 



i=l 
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Therefore for R suitably large enough A maps Vji into itself. In the following we illustrate 
the choice of the constants in detail. First we choose Eq > such that Leq < 1/4 
and Eq < r/2. Then we choose Rq > such that C'{eo) < Ro/2, which means that 
Ro/2 + C'{eo) < Ro- Now for a given arbitrary but fixed R > Rq we choose 6 > such 
that 

CS^R<^ and C{R)6'^ 



where the constants C, C{R) are from inequalities (5.3) and (5.9). With this choice of Eq, 
R and 6 we observe 

C5^R + Eq < r (such that the properties (A1)-(A3) are fulfilled), 

R 

— + C'{eo) < R (such that A is a self mapping), 

^ 1 

C{R)62 + Leq < - (such that A is a 1/2-contraction). 

Therefore we conclude that A has a unique fixed point in V^, which was the remaining 
part to prove the theorem. □ 

Remark 5.2 (A continuation criteria). The question arises on which interval [0, Tmax) the 



mean curvature flow with triple junction (1.1), (1.2) can he extended. A careful revision 



of the above proof shows that 5 in the local existence interval depends on the size of r 
(responsible for the validity of Assumptions (A1)-(A3)) and on Eq. We note that for the 
validity of Assumptions (A1)-(A3) we need that the metric tensor is positive definite and 
in particular that the inverse exists. In Remark \2.1\ we gave a formula for the metric 
tensor and one can see that if the second fundamental form of T\ and terms dirl are 
bounded, we can give a lower bound on the choice of r. If in addition we choose eq small 
enough, this would lead to a lower bound on the existence interval [0,5]. In this way, we 
can achieve existence in any given time interval [0,T] by splitting it into small ones and 
by choosing appropriate reference configurations on each interval, providing the dirl can 
be chosen bounded for all reference configurations on the interval [0,T]. 

We remark that the bound on dirl can be achieved in the following way. If we choose 
the vector tI as a truncation of the unit outer conormal with the help of geodesic lines, we 
can do this in a strip around dTl given by q + ri>l{q), where q G dTl and < r < tq for 
some positive rg. Here we replace r by a cut-off function evaluated at the geodesic distance 
from (9r*. This gives a minimal bound on the diameter of the neighbourhood of the triple 
junction, where rl does not vanish and in this way we can also bound derivatives of the 
form diTl. Possible scenarios for which this cannot be achieved are the following: 

• The area of one hypersurface converges to zero. 

• The triple junction develops during the evolution a self contact. 

A similar continuation criterion in the case of curves has been studied in Mantegazza, 
Novaga and Tortorelli ^25], where the authors consider evolution of planar networks ac- 
cording to curvature flow and conclude existence as long as one of the length of the curves 
tends to zero or a curvature integral blows up at a certain minimal rate. 
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Remark 5.3 (Cluster with boundary contact). We remark that it is also possible to 
consider a configuration where the three hypersurfaces lie inside a fixed bounded region 
Q C M."''^'^ and meet its boundary at a given contact angle, see for example Bronsard and 
Reitich or Garcke, Kohsaka and Sevcovic [T^ for curves in the plane, and Depner 
or Depner and Garcke JE/ for arbitrary dimensions. A natural contact angle achieved by 
the minimization of the weighted area would be 90 degree. If one uses the parametrization 
of or to describe the geometric problem as a system of partial differential equations 
and the ideas from f^, /73| / or from this work, one could derive a local existence result 
also in this situation. 
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